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We discuss spontaneous supersymmetry breaking in the .A" = 1 Wess-Zumino model in two di- 
mensions on the lattice using Wilson fermions and the fermion loop formulation. In that formu- 
lation the fermion sign problem related to the vanishing of the Witten index can be circumvented 
and the model can be simulated very efficiently using the recently introduced open fermion string 
algorithm. We present first results for the supersymmetry breaking phase transition and sketch 
the preliminary determination of a renormalised critical coupling in the continuum limit. 
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1. Introduction 

Supersymmetry as an extension of the space-time symmetries is an interesting concept which 
might well be realised in nature in one form or another. Its presence has several intriguing con- 
sequences, for example, one expects a vanishing ground state energy, if supersymmetry is exact. 
Moreover, the particle spectrum contains mass degenerate bosons and fermions which are related 
by the supersymmetry. However, so far there has not been any experimental sign of such a boson- 
fermion degeneracy. So, if supersymmetry is indeed realised at some high energy scale, it must be 
broken at the low energy scales accessible in today's experiments. The scenario of (spontaneously) 
broken supersymmetry then implies that there is no supersymmetric ground state, the ground state 
energy is not vanishing, and the particle masses need not be degenerate. In this context, it is inter- 
esting to ask how the spontaneous breaking of the supersymmetry is realised. Since spontaneous 
symmetry breaking is an inherently non-perturbative problem one needs non-perturbative methods 
in order to approach it meaningfully. One such method is provided by simulating supersymmetric 
theories on a space-time lattice. However, since the space-time symmetries are explicitly broken 
by the lattice regularisation and are restored only in the continuum limit, also the supersymmetry is 
in general not (or not fully) reaUsed on the lattice. As a consequence, there is a subtle and delicate 
interplay between the various symmetries, and their realisation in the continuum needs to be care- 
fully studied. Here we present preliminary results of such a study for the JV = 1 Wess-Zumino 
model in two dimensions. Using the Wilson fermion discretisation one can formulate the model 
in terms of fermion loops which can be simulated very efficiently using the open fermion string 
algorithm [1]. In addition, the fermion loop formulation provides a way to circumvent the sign 
problem related to the vanishing of the Witten index [2]. 

2. The ,yi^ = 1 Wess-Zumino model 

The ^ = 1 Wess-Zumino model in two dimensions is one of the simplest models which may 
exhibit spontaneous supersymmetry breaking. Its degrees of freedom consist of one real Majorana 
fermion field i/A and one real bosonic field 0, while its dynamics is described by the Lagrangian 
density 



Here, P{^) denotes a generic superpotential, and P' ,P" its first and second derivative with respect 
to (/). In the following we will concentrate on the specific form 



which leads to a vanishing Witten index W = and hence allows for spontaneous supersymmetry 
breaking [3]. The conxsponding action enjoys the following two symmetries. First, there is a single 
supersymmetry given by the transformations 




(2.1) 



(2.2) 



50 = £i//, 8\if = {^(^ -P')e , 5y = 0, 



(2.3) 



and second, there is a discrete Z(2) chiral symmetry given by 




(2.4) 
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where 75 = 03 can be chosen to be the third PauU matrix. 

For the chosen superpotential the Witten index turns out to be zero, as can be seen as follows. 
Integrating out the Majorana fermions yields the (indefinite) Pfaffian PfM of the Majorana Dirac 
operator M. The corresponding partition function with periodic boundary conditions (b.c.) in all 
directions is equivalent to the Witten index, 

y" ^0£'"^b('^)pfMpp(0)ociv, 

where 5'b((/') is the action for the bosonic field. Now, under the Z(2) symmetry (/> — —0 one has 

5b^5b, PfMpp^-PfMpp, 

so for every bosonic field configuration contributing to the partition function, there exists another 
one with exactly the same contribution but opposite sign, hence yielding = 0. This constitutes a 
necessary (but not sufficient) condition for the supersymmetry to be broken spontaneously. In that 
case, one expects a bosonic and fermionic ground state related to each other by the supersymmetry 
transformation. On the other hand, if one chooses thermal b.c. (antiperiodic b.c. for the fermions in 
time direction) the supersymmetry is broken by the finite temperature of the system and one finds 

5b ^Sb, PfMpa^+PfMpa. 

In order to further understand the supersymmetry breaking pattern, i.e. the relation between 
the supersymmetry breaking and the Z(2) symmetry breaking, it is useful to consider the potential 
for the bosonic field, \^P' {^Y = + jg^^'* + const. It simply represents a standard (/)^-theory 

in which, depending on the choice of the bare parameters m and g, the Z(2) symmetry may be 
broken. Indeed, for large values of mjg the Z(2) symmetry is spontaneously broken (in infinite 
volume) and the boson field selects a definite ground state. Denoting with ^ the expectation value 
of the volume averaged boson field, one finds 

= +m/2g ^ Pf Mpp = +Pf Mpa , 
= -mjlg Pf Mpp = -Pf Mpa , 

i.e. in the former case the unique ground state is bosonic, while in the latter it is fermionic. In both 
cases, there is a single, unique ground state tantamount to having unbroken supersymmetry. By 
contrast, for small values of m/g the Z(2) symmetry is unbroken, i.e. one has 0=0 which allows 
both a bosonic and fermionic ground state, tantamount to having broken supersymmetry. Indeed, 
the tunneling between the two equivalent ground states corresponds to the massless Goldstino mode 
which comes along with any spontaneous supersymmetry breaking. 

2.1 Lattice discretisation and fermion loop formulation 

In order to put the Wess-Zumino model on the lattice we follow the approach of Golterman 
and Petcher [4] where it is shown that using the same lattice derivative for the bosons as for the 
fermions (and renormalising the mass parameter m accordingly), the supersymmetry is guaranteed 
to be restored in the continuum limit. 
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Using the Wilson lattice discretisation for the fermion fields yields the fermion Lagrangian 
density 

where <^ is a real, 2-component Grassmann field, ^ = — is the charge conjugation matrix and 
5*, 5 are the backward and forward lattice derivatives, respectively. However, while the Wilson 
term d*d avoids fermion doubling, it spoils the discrete chiral symmetry of the fermion action as 
well as the Z(2) symmetry — )• —0 of the boson action^ 

Another problem for simulating the model on the lattice is the fact that the Pfaffian is indefinite. 
As discussed above, this is due to the vanishing of the Witten index and constitutes a generic 
problem for any numerical Monte Carlo investigation of spontaneous supersymmetry breaking, 
independent of the chosen discretisation. The problem stems from the fact that the indefinite Paffian 
can not be simulated directly by standard Monte Carlo methods. Instead one uses the effective 
action 

5eff(0)=5b(0)-ln|PfM(0)| 

for the boson field and takes the sign of the Pfaffian into account by reweighting. This approach 
in general leads to severe sign problems [5,6]. 

It turns out that the sign problem can be circumvented for Wilson fermions by using an ex- 
act reformulation of the lattice model in terms of fermion loops [2] and simulating fluctuating 
fermionic boundary conditions [1]. In the loop formulation one expands the Boltzmann factor of 
the fermion action, effectively constructing a hopping expansion. When one subsequently per- 
forms the integration over the fermion fields, the nil-potency of the Grassmann elements ensures 
that only closed, non-oriented and non-intersecting fermion loops survive. The partition function 
then becomes a sum over all self-avoiding fermion loop configurations I, 

{£'}Gi? 

where G)[£, 0] denotes the weight for a given loop configuration I, and denotes the equivalence 
class of loop configurations with an even or odd number of loops winding around the lattice in the 
spatial and temporal direction, respectively. Z_sf represents a system with unspecified fermionic b.c. 
[7], while the system with periodic b.c. for the fermion, i.e. the Witten index, can be constructed 
by forming 

W EE Zpp = z^o„ - z^jo - Z^„| - z^„ , 
or the system at finite temperature by forming 

Note that the weight (O does not necessarily need to be positive definite in each of the sectors, but 
in practice it turns out that it is the case as long as one stays close enough to the continuum limit. 

As described in [1] the system can most efficiently be simulated by introducing an open 
fermion string corresponding to the insertion of a Majorana fermion pair. By letting the ends 

'in principle this complication can be avoided by using a lattice discretisation which respects the discrete Z(2) 
chiral symmetry, e.g. the SLAC derivative [5]. 
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Figure 1: Monte Carlo time history of at mjg = A,ag = 0.125 (left plot) and m/g = 0.16, ag = 0.03125 
(right plot). 

of the string move around the lattice by a standard Metropolis update procedure, one samples the 
fermion 2-point function as well as the relative weights between Z , Z j^^g , Z^^j and Z^j , . Finally, 
the bosonic fields are simulated by standard Monte Carlo methods. 



3. Results 

By looking at the behaviour of for large and small values of m/g one can check for the Z(2) 
symmetry breaking^. The plots in Fig. 1 show the Monte Carlo time history of atm/g = 4,ag = 
0.125 (left plot) and m/g = 0.16, ag = 0.03125 (right plot). For m/g = 4 the system is in the Z(2) 
broken phase with {\<p\) ~ 2. Since the system is at finite volume it still tunnels between the two 
vacua with (|^|) ~ ±m/2g, but the tunnelling will be suppressed in the limit L — )• oo or m/g — > oo. 
For m/g = 0.16 on the other hand, the system is in the Z(2) symmetric phase with (0) ~ 0. 

It is now interesting to see how the partition functions Zjfgg,Z^^g,Zjfg^ ,Zjc^^, or Zpp and Zpa, 
behave in the two situations. Fig. 2 shows the probability distributions of the various sectors as a 
function of 0, again for m/g = 4,ag = 0.125 (left plot) and m/g = 0.16, ag = 0.03125 (right plot). 
In the first situation where the Z(2) symmetry is broken, one finds 

((/») ~ -2 : Zoo ^ Zio ~ Zoi ~ Zii Zpp ~ -Zpa , 

((^) ~ +2 : Zoo ^ 1, Zio ~ Zoi ~ Zii ~ Zpp ~ +Zpa , 

so (0) ~ — 2 corresponds to the fermionic ground state while (0) ~ +2 coiTcsponds to the bosonic 
one. In either case, a unique ground state is chosen by the system (up to finite volume tunneling) 
and hence supersymmetry is unbroken (at least in the thermodynamic and continuum limit). 
In the second situation, where the Z(2)-symmetry is unbroken, one finds 

(0)~O: Zoo ~Zio + Zoi +Zii Zpp^^O, 

so the bosonic and fermionic ground states occur with equal probability (thereby cancelling their 
contribution in Zpp = W) and hence supersymmetry is spontaneously broken. 

As a next step one can now determine the value of m/g at which the transition from the 

^Note that (j> is not a true order parameter for the Z(2) symmetry since the symmetry is explicitly broken by the 
lattice discretisation. 
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Figure 2: Probability distribution of the partition functions Z^^^ , Z^^, , Z^jf^j I'^^xx ^ ^ ^^"^ = 4, flg = 
0.125 (left plot) and mjg = 0.16, ag = 0.03125 (right plot) as a function of ^. 

Z(2) symmetric and supersymmetry broken phase to the Z(2) broken and supersymmetric phase 
occurs. This is most easily done by scanning am at fixed lattice spacing ag for various lattice 
extents L. The critical value anic where the phase transition occurs determines the dimensionless 
critical coupling m^jg at the given lattice spacing. The procedure is illustrated in Fig. 3 where 
the left plot shows the ratio Zpp/Z serving as a (pseudo-)order parameter for the supersymmetry 
breaking phase transition, while the right plot shows (i'(^)pa = (sign^)pa as a (pseudo-)order pa- 
rameter for the Z(2) symmetry breaking^, as a function of the bare mass at fixed lattice spacing 
ag = 0.03125 for various lattice extents. We note that the behaviour of seems to suggest a 
second order phase transition. Setting up the model with the Wilson derivative for bosons and 
fermions yields a supersymmetric continuum limit [4]. Since the model is superrenomalisable it is 




Figure 3: Ratio Zpp/Z (left plot) and {s,fi) ~ (sign^) serving as a (pseudo-)order parameter for the super- 
symmetry breaking phase transition and the Z(2)-symmetry breaking phase transition, respectively, as a 
function of am at fixed lattice spacing ag = 0.03125 for various lattice extents. 



^Note that since Zpp ~ in the supersymmetry broken phase, expectation values need to be calculated in the thermal 
ensemble in order to be under good numerical control. 
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sufficient to tune only the mass parameter m 
in order to obtain a renormalised theory. Us- 
ing the relation = m| + 2^^/71 In m| be- 
tween the bare mass m and the renormalised 
one m^, and setting the scale, i.e. the lat- 
tice spacing a, by the dimensionful coupling 
g = g/a, one can determine the continuum 
hmit of a dimensionless critical coupling for 
the supersymmetry breaking transition with 



/c. 



lim — 

g^O Mr 



The procedure is illustrated in Fig. 4 and it 
will be interesting to see how this result com- 
pares to previous determinations [5,8]. 




0.15 



Figure 4: Continuum limit of the renormalised di- 
mensionless critical coupling /ait for the supersym- 
metry and Z(2)-symmetry breaking phase transition. 



4. Outlook 

There are a several obvious ways to proceed from this first, preliminary investigation. Firstly, 
one can use Wai^d identities as (pseudo-)order parameters to determine the phase transition point. 
Secondly, one can determine the boson and fermion mass spectra. The latter is particularly simple 
in the fermion loop formulation. Thirdly, one can try to perform non-perturbative renormalisation 
using the boson or fermion masses. Finally, it would also be interesting to implement the loop 
formulation of the model with a domain wall or overlap type fermion discretisation for which the 
discrete Z(2) chiral symmetry remains exact at finite lattice spacing. 
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